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ABSTRACT 

We implement the theory of resonant scattering in the context of Cosmic Microwave 
Background (CMB) polarisation anisotropies. We compute the changes in the E-mode 
polarisation (EE) and Temperature E-mode (TE) CMB power spectra introduced by 
the scattering on a resonant transition with a given optical depth tx and polarisation 
coefficient E\. The latter parameter, accounting for how anisotropic the scattering 
is, depends on the exchange of angular momentum in the transition and enables, 
a priori, to distinguish among different possible resonances. We use this formalism 
in two different scenarios: cosmological recombination and cosmological reionisation. 
In the context of cosmological recombination, we compute predictions in frequency 
and multipole space for the change in the TE and EE power spectra introduced by 
scattering on the H a and P a lines of Hydrogen. This constitutes a fundamental test 
for the standard model of recombination, and the sensitivity it requires is comparable 
to that needed in measuring the primordial CMB B-mode polarisation component. In 
the context of reionisation, we study the scattering off metals and ions produced by the 
first stars, and find that polarisation anisotropies, apart from providing a consistency 
test for intensity measurements, give some insight on how reionisation evolved: since 
they keep memory of how anisotropic the line scattering is, they should be able to 
discern the OI 63.2/im transition from other possible transitions associated to OIII, 
Nil, NIII, etc. The amplitude of these signals are, however, between 10 to 100 times 
below the (already challenging) level of CMB B-mode polarisation anisotropies. 

Key words: cosmic microwave background - galaxies: clusters: general - diffuse 
radiation - intergalactic medium 



1 INTRODUCTION 

The cosmic microwave background (CMB) constitutes one 
of the most powerful tools in cosmology today. Observations 
of the angular power spectrum of the CMB intensity and 
polarisation anisotropies are providing constraints at the 
level of few percent on the values of the cosmological 
parameters describing our Universe with in the standard 
cosmological scenario, (e.g. ISpergel et ail (120061 )1. More- 
over, because these anisotropies are intricately linked to the 
initial conditions in our Universe, they can also be used as 
a probe of inflation and fundamental physics. Microwave 
intensity (or temperature) anisotropies have recently been 
measured in the whole sky from large to medium (~ 0.2°) 
angular scales by the Wilkinson Microwave Anisotropy 



Probe, (WMA10). WMAP measurements have been com- 
plemented by othe r ground-based high-angular resolution 
CMB experiments (Ipickinson et~all 12004 iKuo et all 120041 : 
iReadhead et al.ll2004n " with partial sky coverage. 



On the other hand, the study of the CMB polarisa- 
tion anisotropies is ju st starting: experiments like CAPM AP 
jBarkats et all 120051 ) and CBI l|Readhead et all 12004 ) an- 
nounced the first detection of CMB polarisation anisotropies 
in subdegree angular scales, whereas WMAP's team man- 
aged to separate large-scale CMB polarisation anisotropies 
from the dominant contamining signals produced in our 
Galaxy. Indeed, when trying to achieve better constraints 
on the theoretical models and to test the inflationary pre- 
dictions, polarisation measurements are known to play an 
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important role in breaki ng some of the de generacies present 
in the parameters (e.g. IPage et al.l 120061 ). Furthermore, a 
particular type of CMB polarisation anisotropies having 
particular symmetry properties (the so-called B modes of 
polarisation) are known to be closely related to tensorial 
modes in the metric perturbations (that is, gravitational 
waves) generated during the inflationary epoch. Therefore, 
it is easy to understand that the search for these polarisation 
anisotropies provides a unique test for both the cosmological 
model and fundamental physics, since to date no detection 
of gravitational waves has been reported. 

For this reason, in the last few years an important ef- 
fort has been put into measuring the polarisation angular 
power spectrum of the CMB, and consequently there is al- 
ready a considerable number of projects (Clover, QUIET, 
CMBPol) pursuing the detection of the primordial B-mode 
of the CMB polarisation anisotropies. The predicted sensi- 
tivities of these experiments (10 -2 to 10 -4 (fiK) 2 ) will open 
the possibility to investigate additional interesting effects 
whose amplitudes were regarded as unreachable. 

In particular, the resonant scattering on Hydrogen 
Balmer and Paschen lines at the epoch of recombination 
[z ~ 1050) is known to produce frequency-dependent tem- 
peratur e anisotropies in the sub-deg ree angular scales of the 
CMB, l|Rubino-Martm et all [2005T ). These same lines are 
known to give ris e to spectral distortions in the CMB black 
body spectrum l|Rubino-Martm et alj 120061 ; IChluba et al.l 
2006), at the relative level of ~ 10 _s , and in general of 
the order of ~ 10~ 7 for higher series. Here we analyse the 
changes that these transitions introduce in the polarisa- 
tion angular spectrum, and make definite predictions for 
the standard model of recombination. The amplitudes of 
the changes in the polarisation power spectrum are found 
to be comparable to the typical amplitudes of the B modes 
for values of the tensor-to-scalar ratio r ~ 10 -2 — 10~ 3 . 
These signals should constitute a new probe for the end of 
recombination process, and may provide information about 
cosmological parameters and their evolution with time. 

Another physical process that upcoming CMB experi- 
ments may be sensitive to is the resonant scattering off met- 
als and molecules at high redsfhift. This idea is not new, and 
several authors have studied the effect of molecules (LiH, 
DH + ) on the CMB radiation a t high redshift, both from 
the theoretical jDubrovichl 19771. 1993) or the o bservational 
point of view, I de Bernardis et al ] Tl993h . In iBasu et al.l 
(|2004 l (hereafter BHMS), a detailed theoretical study of the 
effect of resonant scattering of CMB photons off metals and 
ions was presented. In that work, it was shown that different 
observing frequencies would probe the resonant scattering at 
different redshifts. For this reason, the comparison of CMB 
maps obtained at different frequencies should provide infor- 
mation of the metal enrichment history of the Universe. This 
is essentially a consequence of the frequency dependent op- 
tical depth (r„) that these resonant transitions give rise to. 
In their analysis, BHMS studied fine structure transitions 
of atoms and ions (CI, CII, Nil, NIII, OI, OIII, etc), and 
found that the scattering process should be dominant over 
collisional e xcitation. Further, they exploited the fact that, 
as shown in IHernandez-Monteagudo fc Sunvae 3 l|2005l ). the 
changes in the CMB temperature angular power spectrum 
are linear in t v , and proposed a method to extract the metal 
induced signal in upcoming CMB data sets. 



IHernandez-Monteagudo et alj l|2006T l addressed the re- 
quirements that future CMB data must fulfill in order to be 
sensitive to the metal-induced signal. Among all of them, the 
most demanding was the cross-channel calibration, whose 
relative uncertainty had to be below t„, i.e., better than 
~ 10~ 4 . Point spread function reconstruction errors were 
found not to be very critical. The observing channels should 
be placed in an optimal frequency range of [100, 350] GHz, 
enabling an accur ate subtraction of the ther mal Sunyaev- 
Zel'dovich effect, (|Sunvaev fc Zel'dovich|[i980h . 

In this paper, we study the changes that resonant scat- 
tering off metals introduce in the CMB polarisation angular 
power spectrum. The amplitude of these changes is found to 
be extremely small, still a factor 10 - 100 below the nominal 
sensitivity of proposed B-mode experiments. However, they 
would provide a parallel test for the enrichment history of 
the Universe, and help to distinguish the OI among other 
likely resonant species, like Nil, NIII, CI, CII or even OIII. 

This study should be understood in the context of pro- 
posed CMB polarisation measurements, that attempt to 
measure primordial B-modes of typical amplitudes 10 -3 
(fiK) 2 . Distinguishing signals of such small amplitudes de- 
mands stringent requirements not only on the sensitivity 
side, but also (and perhaps more importantly) on the sys- 
tematic/foreground side: contaminants are now known to 
be orders of magnitude larger that the targeted B-modes, 
and their exquisite control and removal will be absolutely 
demandatory. The signals we are addressing in this paper 
have comparable amplitudes to the intrinsic B-modes: if the 
latter are finally measured, then these signals, containing a 
wealth of information about the recombination and reioni- 
sation histories of our Universe, may be accessible too. 

The structure of the paper is as follows. In section 2 
we implement the theory of resonant scattering in the con- 
text of CMB angular anisotropies. In Section 3 we apply 
our formalism to the polarisation changes introduced by the 
resonant scattering on H a and P a lines of Hydrogen during 
cosmological recombination. In Section 4 we address the res- 
onant scattering on metals and ions during the end of the 
Dark Ages and the beginning of reionisation. We discuss our 
results and conclude in Section 5. 



2 RESONANT SCATTERING AND CMB 
POLARISATION ANISOTROPIES 

2.1 Background 

We start by reviewing well known results on resonance radi- 
ation and CMB polarisation. These are combined in the next 
subsection in order to compute the changes in the CMB po- 
larisation power spectrum for different types of transitions. 
The reader familiar with the formalism introduced here may 
therefore prefer to skip this subsection and continue in the 
next one. 

An analysis of the polarisation state of resonance radi- 
ation com i ng ou t of a generic atomic system is presented in 
iHamiltonl l| 19471 ) and IChandra sckhar (1960). Their main re- 
sult is that the resonant scattering is very close to Rayleigh 
scattering, although in the general case it gives rise to a 
smaller degree of polarisation. Indeed, the resonant scatter- 
ing can be decomposed into a Rayleigh-like scattering term 
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(i.e., with its same angular dependence) and an isotropic 
scattering term. If 7m and 7^ denote the incident light in- 
tensities parallel and perpendicular to the scattering plane 
in a solid angle element around direction n', then the scat- 
tered counterparts (7y and I±) along direction n are given 
by 



7x 



dfi(n) = tx 



F 3 



E 2 — 

Sir 



1 1 
1 1 



7 



dn(n'). 



(1) 



8 and rx are the scattering angle and optical depth, re- 
spectively. The coefficients E\ and E2 provide the relative 
weight of the Rayleigh-like and isotropic scatterings, and 
verify the condition E\ + E2 = I , whereas the other co- 
efficients (proportional to 1/471") ensure conservation of the 
scattered total intensity. The values of E\ and E2 depend on 
the initial value and change of the quantum number describ- 
ing the angular m omenta involved in the resonant transition, 
|Hamiltorjll947T ). 

We next write the radiative transfer equation for this 
scattering in terms of the Stokes parameters I (= 7m + 7_l) 
and Q (= 7ii — 7x). This scattering produces no circular 
polarisation, so the Stokes parameter V will be ignored. We 
study the radiation coming from an arbitrary point n on the 
celestial sphere. Our axes used to define the parallel (||) and 
perpendicular (_L) directions are determined by the merid- 
ian and azymuthal dir ections at n, respectively. As shown in 
IChandrasekhar] l|l960l) . in order to compute the parameters 
I(n) and Q(n) in this reference frame, it is necessary to in- 
tegrate in all possible directions n' of incoming radiation. In 
this step, one must account for the different reference frames 
being involved. Since our scattering matrices (eq.{T}) are de- 
fined with respect to the plane of scattering, it is necessary 
to make two transformations of the intensities: the first one 
projects intensities at n' into the axes parallel and normal 
to the scattering plane, whereas the second transforms the 
scattered intensities into our reference frame at n. Under 
an infinitesimally small optical depth due to the resonant 
transition (Arx), the radiative transfer equation yields the 
following changes in the I and Q parameters for the resonant 
scattering: 



(n) = Arx - * (n) + 



7 - Et P2M^(h + Q2 + Qo) 
Ex (l - P 2 (/i))§(7 2 + Q2 + Qo) 



(2) 



In this equation, the parameters I and Q have been expanded 
on a Legendre polynomial basis (X(ii) = ^ ; (— i)'(2/ + 
l)Pi (fi)Xi, with X = I,Q), [i = cos 8 and 9 is the polar angle 
at n. Both I and Q are assumed to be azymuthally symmet- 
ric. This expression meets the standard Rayleigh scattering 
limit if 75 1 — > 1, and shows how the intensity quadrupole 
alone (72) acts as a source for linear polarisation. At the 
same time, it is clear that different resonant transitions with 
different Ei coefficients give rise to different levels of polar- 
isation. We shall show below that this fact may have its 
relevance, since polarisation measurements may distinguish 
among different candidate resonant species otherwise undis- 
cernible with intensity measurements only: the polarisation 



term II = (I2 + Q2 + Qo) is typically a few percent of 7o, 
and this fact largely hides the value of Ei from intensity 

observations. 

In a cosmological context, we follow iMa fc Bertschingeil 

(1995) to describe the evolution of the perturbations of the 
photon field. They work in Fourier space, and therefore the 
photon distribution function / 7 (k, q, 77) is a priori dependent 
on comoving wavemodes k, photon momenta q and confor- 
mal time 77. However, in cosmological perturbation theory, 
due to a symmetry present in the evolution equations, it is 
customary to work with the simplifying assumption that / 7 
depends on the direction of the photon momentum q/q = n 
via the dot product nk = n-k/fc only, (k being the modulus 
of k and q the modulus of q). This permits to expand the 
photon distribution function for the polarisation state i as 



/ 7 (k,q,r?) = fo(q,v) 



fo(q,v) 



'(2 + l)fl(k.n)¥f(k,g,ij) 



(3) 



with /o(q, 77) = l//i 3 (exp (qc/kBTcMB[rf\) — 1) the unper- 
turbed Planck distribution function for the polarisation 
state i. "J; is the vk, q and rj dependent perturbation of 
the distribution function for the same polarisation state i. 
This expansion can be applied to the momentum-averaged 
energy distribution function for both polarisation states, 



ELi / d ii 2 fo(q>v)<i 

and for the difference of the two polarisation components 



(4) 



G 7 ,;(M) = 



/ dqq 2 q(f^q, V )<fl - fS(q, V )<f 



(•>) 



These two /-dependent quantities correspond to the Xi 
terms in the expansion A(k,q, 77) = ^ ; (— i)'(2Z + l)P;(k ■ 
n)Xi(k, q, rj), with X = F 7 ,G 7 . If now a reference system 
for q is chosen so that the z axis is parallel to k, then it is 
clear that the linear polarisation will be in the direction of 
changing 0, i.e., parallel to d/88, due to the azymuthal sym- 
metry of / 7 (k, q, k • n, 77) for fixed k and 77. Therefore, in the 
orthonormal basis {d/88, d/d(j>} the momentum-integrated 
quantity G 7 corresponds to the Q Stokes parameter, whereas 
U vanishes. IZaldarriaga fc Hararil l|l995l ) formally integrated 
the evolution equations for At = and Ap = G 7 /4 

for every Fourier mode k, (note that we have switched into 
temperature units). However, since the Q and U Stokes pa- 
rameters are dependent on the orientation of the reference 
frame where they are defined, when integrating in all Fourier 
modes to obtain At, Ap along a given line of sight n the 
following problem arises: for every new k a new (rotated) 
reference system must b e adopted, and Q and U m ust be 
transformed accordingly. Zal darriaga fc Seliakl (119971 ) solved 
this problem by introducing the combinations Q(n) ± ill (n) 
that can be easily transformed into rotationally invariant 
quantities E(n) and 73(n). After projecting the three quan- 
tities Ar(n), E(h) and B(h) on a spherical harmonic basis, 



A t ,e,b = ^2 a l,m' BY l,m(n), 



(6) 
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they found analytical expressions for the angular power 
spectra C ; T 



yT,E,B I T.E.B I T,E.B\*\ 

Km ) r- 



c 



(4tt) 2 / dk k 2 P,p(k)\Af' E '(Mo)!' 



(7) 



P<p(k) is the initial power spectrum of scalar perturbations, 



which are found not to give rise to any B mode, (C; 
The transfer functions 



'(Mo) 



drj S T ,E(k,7i)f l ' ji(k[r) - 77]) 



0). 

(8) 



are characterised by the source functions 

S T (k,rj) = e~ T (i + a) + a(a t ,o + a J + A( a+ ^ 



k \ 4 



for the temperature and 

3An 
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777T + A THT + A 77T 
4fc 2 / 4fc 2 4fc 2 



(9) 



S B (M) = 



4(fc[r? - V]) 2 



(10) 



for the E-mode polarisation field. /; 



1 VZ and /, 



y/ (I + 2)!/(Z — 2)!. In these equations, ji(x) are spherical 
Bessel functions of order I, 770 is the conformal time at 
the present epoch, and a = (£ + Qr/)/2k 2 is a function 
of the metric scalar perturbations h and £ in the syn- 
chronous gauge, used throughout this paper. The polarisa- 
tion term II, as already defined above, is the sum of the 
temperature quadrupole and Q monopole and quadrupole, 
i.e. II = At,2 + Ap.2 + Ap,o, and A = re _T is the visi- 
bility function, f is the opacity, which, in case of Thomson 
scattering, reads 



(11) 



with a the scale factor, ot the Thomson cross section and 
n e the comoving free electron number density. The optical 
depth is defined as 



t t (v) 



f T (v')dv'- 



(12) 



For these quantities, the subscript T denotes "Thomson 
scattering". Throughout this paper, dots denote derivatives 
with respect to conformal time rj. 

2.2 The Modified Visibility Function 

We combine the previous results in order to compute the ef- 
fect of resonant scattering on CMB anisotropies. In equation 
we have seen that a resonant transition of optical depth 
Atx erases both I and Q parameters by an amount Ar rs , 
and generates new intensity by an amount proportional to 
Atx as well. However, the new polarisation generated in 
the resonant scattering is not proportional to Atx, but to 
Ei Atx- This fact requires the following modifications of 
eqs. (|9HUp by means of the new visibility function A: 



S T (k,r)) = e T (i + a)+A A T ,o + ^ +A [a + 
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Figure 1. Function H displaying the ratio of the drag force 
exerted by resonant scattering over the Thomson drag, for two 
different resonant transitions: Ly-ct (solid line), and Lil 2S^ 2P 
[6708 A] (dashed line). At all redshifts, S is much smaller than 
one, making the drag force induced by resonant species negligible 
if compared to the Thomson drag. 



q v 3An 

In this case, the total opacity and optical depth read 
f(r?) = T T {n) + Tx{n), 



-to) = 



(tt + T X )dr/'. 



(14) 

(15) 
(16) 



The visibility function A remains unchanged in terms of r 
and f, while the new visibility function A is defined as 



A = (tt + Eiix) 



(17) 



Note that this new visibility function is coupled to the po- 
larisation terms exclusively. The generation of anisotropies 
is ruled by the opacity term, which does contain E\ , whereas 
the blurring of existing anisotropies is caused by the expo- 
nential of r, which does not depend on E\. Hence, we expect 
that different resonant transitions with different E\ coeffi- 
cients should be discernible in those angular scales (I ranges) 
where the generation of new polarisation anisotropies dom- 
inates or is comparable to the blurring of old anisotropies. 
I.e., we shall recover no information on E\ in the high I limit 
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2.3 The Drag Force 

The Thomson scattering influences the evolution of baryons, 
since they exchange momentum with the CMB photons. 
This interaction is accounted for by adding a momentum 
conservation term to the equ ation for the evolution of t he 
baryon peculiar velocity, (e.g.. lMa fc Bertschingerl l|l995h V 
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-Vb 



ic^kSb + 1 
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(18) 



3p b v 4 

ff, denotes baryon peculiar velocity, c 3 is the sound speed 
in the baryon photon fluid before recombination, 5b is the 
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comoving baryon density contrast, and PjtM are the photon 
(baryon) background energy densities. Finally, F 7] i is the 
I — 1 moment in the expansion of the momentum averaged 
photon energy distribution function, (see eq.Q). When con- 
sidering the resonant scattering, a similar term should be 
included in this equation. There is, however, a critical dif- 
ference: Thomson scattering is frequency independent, and 
couples all photons with baryons, as opposed to resonant 
scattering, that couples only t hose p hotons within the ther- 
mal width of the line. lLoebl |200lf ) uses arguments based 
upon the characteristic time of the drag force to conclude 
that it does not play any role in the case of Lithium scat- 
tering. Here we explicitely compute the drag force, for both 
Thomson scattering and resonant scattering, and confirm 
that indeed the latter can be ignored in practically all cases. 

For Thomson scattering, the drag force (defined here as 
rate of momcmtum exchange per unit time and unit volume) 
reads 

F 7 T (k )J7 ) = (-47ri) J dqq 2 J2fo(q,n)qViQt,q,n) x 

i — 1 

ca T n e (rj) (19) 

but for resonant scattering on a transition X, the drag 
force must acccount for limited range of frequencies that do 
interact with the baryons. This is introduced via by impos- 
ing a resonance condition in the momentum integration of 
the atom distribution function: 



2 

F 7 X (k,77) = (-4tu) / dqq 2 J2fo(q, vWifcw) X 

=i 



' ' t y / dpf x (p,v)- 



(20) 



The symbol E denotes the subset of all possible atom mo- 
menta p verifying the resonance condition. fx(p,rj) ' s t ne 
distribution function of the atom responsible for the tran- 
sition X. For simplicity, we shall assume that all baryons 
are in thermal equilibrium with the CMB radiation field 
at all redshifts, so that the CMB temperature will charac- 
terise their Maxwell-Boltzmann distribution function. The 
resonance condition 



_ ( P n 

h \mxc 



qc 



- 1 = v x 



(21) 



can be rewritten, after accounting for the natural width of 
the line, as 



p ■ n £ 



l ( mXC ( 1 ~? ± lf) 



max [ mxc | 1 ± — — 

v 4-7T 



(22) 



Ax is the Einstein coefficient for spontaneous emission of 
transition X, and the ratio Ax /27r is its natural width, mx 
and vx denote the mass of the resonant atomic species and 
the resonant frequency, respectively, p is the cosine of the 
angle between p and n, (p = p/p • n). Since we are con- 
sidering all photon frequencies, all atoms within the ther- 
mal width are counted and all of them contribute to the 
drag force, so the integral on the atom distribution function 



yields nxirj), the atom number density, (see Appendix A for 
details). Therefore: 



F 7 X ~ (-Am)(ca x n x ) ^q 2 x (/o*i) 9 



qx 



(23) 



I.e., momentum of amplitude qx = hv x /c is transferred 
to the atom fluid at a rate per unit volume c oxnx 
in a momentum-space shell of radius qx and thickness 
hvx /(2nc). The photon distribution momentum /o^i corre- 
sponds to the sum over the two polarisation states, /o^i = 
X/i-i /o^i- The ratio of the Thomson induced drag force 
and the one generated by resonant scattering reads 



oxnx 



^7<7x {fo(qx,ri)^i(k,q x ,ri)) qx 



dq q /o(g, J?)*i(k,q,7/) q 



(24) 



If we now consider that, at least ab initio, the / = 1 momen- 
tum of the photon distribution function fo(q, 7?)\&i(q, q, n) is 
caused by gravitational sources, so that d^i/dq ~ 0, then 
$i is cancelled in eq. (|24[l . yielding 
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The ratio S versus redshift is shown in Figure (JTJ for two 
resonant transitions: Ly-a (solid line) and Lil 2S^ 2P [6708 
A] (dashed line). There are two factors that make the drag 
force exerted by resonant transitions negligible when com- 
pared to Thomson drag: i) Only photons within the line 
width exchange momenta with the atoms, and ii) the line 
is, for most redshifts, far from the peak of the planckian 
photon distribution function. Therefore, for the moderate 
to small values of tx that we will be handling in this paper, 
the drag force can be completely ignored in our calculations. 



2.4 Moderate and Low tx Examples 

In this subsection we compute the changes in the intensity 
and polarisation CMB angular power spectra introduced by 
a resonant transition in two different limits: moderate opti- 
cal depth (tx = 2), and small optical depth (tx = 5x 10 -4 ). 
As shown in BHMS, in the limit of low tx there are two com- 
peting terms: a generation term, present at scales larger than 
the projected horizon at the epoch of scattering, and a blur- 
ring term, accounting for the slight decrease in the amount 
of previous-to-the-scattering temperature anisotropics. The 
latter is dominant at small angular scales, and veryfies 
8Cf T ~ — 2tx Cj T . The former term is caused by the 
growth of peculiar velocities with cosmic time, and is respon- 
sible for positive 5Ci 's below a multipole lx ~ (vo~nx)/rio, 
()Hernandez-Monteagudo &: Sunvaevl 120051 ). In the limit of 
high tx, however, the generatio n of new anisotropics dom i- 
nate in the entire angular range, (|Zaldarriaga fc Loebll2002l ) . 

We consider a resonant transition at an effective scat- 
tering redshift close to zx — 500, or conformal dis tance 
rjx — 450 Mpc. This is close to the epoch at which lLoebl 
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Figure 2. TT, EE and TE angular power spectra for three different scenarios: standard ACDM cosmology (thick solid line), standard 
ACDM cosmology plus a resonant transition placed at zx — 500 and (tx, E\) = (2, 1/3) (dashed line), and standard ACDM cosmology 
plus a resonant transition placed at zx — 500 and (tx,Ei) = (2, 1) (dotted line). 
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Figure 3. Relative increment of TT, EE and TE angular power spectra with respect to the standard ACDM scenario due to the presence 
of a resonant line placed at zx — 500 with (tx , E\) = (5 X 10 -4 , 1/3) (dashed line) and (tjc,-Ei) = (5 X 10 -4 , 1) (dotted line). Note that 
for panels (a) and (b) we are plotting absolute values. Since the blurring of original anisotropies is independent of Ex, both lines converge 
to 2tx in the high-/ range for the EE plot as well. Note that due to the change of sign of Cf B , in panel (c) we prefer to normalise by 



l|200ll ): [Zaidarriaga fc Loebl (|2002h placed the resonant scat- 
tering between the CMB and the neutral Lithium via the Lil 
2S^ 2P [6708 A] resonance. Although it has been shown 
that at such redshifts Lithium should b e ionised by Ly-a 
photo ns generated during recombination l)Switzer fc Hiratal 
2005), we choose a similar scattering redshift for compar- 
ison purposes. We use a modified version of CMBFAST 
l|Seliak fc Zaldarriaga|[l996l ) where a new visibility function 
(modified as outlined in subsection l|2.2)l l has been intro- 
duced. Motivated by our results above, we also disconnect 
the drag force exerted by the atoms. We choose a gaussian 
profile for the line opacity, 



tx(v) 



CXD ^-^x) 



Tx- 



(26) 



The width of the line is taken to be one percent of r\x , but 
it would ultimately depend on the spectral resolution of the 
detector. Figure ([2]) shows the temperature (TT), E-mode 
polarisation (EE) and cross (TE) power spectra for three 
scenarios: the concordance ACDM model (thick solid lines), 
used as a reference model, and this reference model as seen 
after introducing a species giving rise to a resonant scatter- 
ing of parameters [tx, Ei]) = [2, 1/3] (dashed lines) and [tx, 
Ei]) = [2,1], (dotted lines). Panel (a) shows that different 
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Table 1. Possible H c 
transitions ( a ) 



(n 



3) and P a (n ■■ 



4) 



Wavelength 




i 


5 


Ei 




(A) 


(s-l) 


{Lj) 


(Lj) 




H 



6562.7096 
6562.7247 
6562.7517 
6562.7714 
6562.8516 
6562.8672 
6562.9093 



5.388e+07 
2.245e+07 
2.104e+06 
2.245e+07 
6.465e+07 
1.078c+07 
4.209e+06 



Pi/2 

Sl/2 
■Pi/ 2 
Sl/2 
P 3/2 
P 3 /2 
PS/2 



D 



3/2 
M 3/2 



Si 



/2 
-1/2 
5/2 
3/2 



Sl/2 



0.5 
0.5 




0.28 
0.32 





18750.684 
18750.720 
18750.829 
18750.881 
18751.011 
18751.011 
18751.064 
18751.065 
18751.111 
18751.138 
18751.191 
18751.210 
18751.225 



5.864e- 
3.065e- 
6.117e- 
3.065e- 
7.037e- 
1.287e- 
3.475e- 
1.173e- 
1.379e- 
9.193e- 
3.128e- 
1.223e- 



06 
06 
-05 
06 
06 
-07 
-04 
-06 
-07 
-05 
-05 
06 



Sl/2 
Pl/2 
Sl/2 
P3/2 



D 



3/2 
P3/2 

Sl/2 
Pl/2 



D 



5/2 



3/2 
'3/2 



"5/2 
P 3/2 



3/2 



D 6 



D 



3.475e+05 D 



/2 
/2 
5/2 
^3/2 
3/2 



D 3/2 
F7/2 
F5/2 
P3/2 
Sl/2 
Pl/2 



0.50 
0.50 




0.28 
0.28 
0.32 
0.32 
0.21 
0.37 
0.02 







(a) Taken from NIST URL site: 

http : //physics .nist .gov/PhysRef Data/ASD/lines_f orm. html 



values of Ei are not discernible in the TT power spectrum, 
since the contribution of the polarisation terms is of the or- 
der of a few percent: in the plot these differences are visible 
to the eye at I < 20, but these are the scales more affected 
by Cosmic Variance. Nevertheless, the EE power spectrum 
shown in panel (b) provides very different estimates for dif- 
ferent values of Ei , making both scenarios easily discernible 
for this large value of tx- Since new anisotropics are gener- 
ated in an older universe subtending a larger angle in the 
sky, the angular pattern of the power spectra shiftes to lower 
multipoles, both in the TT and the EE cases. The polari- 
sation source term SE{k,rj) is larger at z ~ 500 than at 
decoupling (see eq. (|14p l. and this makes the EE power spec- 
trum of the [tx, Ei]) = [2,1] of higher amplitude (dotted 
line above solid lines) . Note that the cross-power spectra also 
show significant differences for the two Ei values considered, 
(Fig.©)). 

The low tx case, displayed in Fig. ([3]), also provides 
significant differences between the Ei = 1 and the Ei = 1/3 
models in the TE and EE power spectra. These differences 
vanish in the high / limit, where the blurring term takes 
over and SCi — > —Itx Ci for all values of E\. fn case of 
resonant scattering on metals produced at late epochs (end 
of the Dark Ages and reionisation) , the generation of new 
anisotropies arise at very large angular scales, and hence 
distinguishing different values of Ei will require large sky 
coverage. 



0.20 L 



1200 
Redshift 



Figure 4. Polarisation coefficients Ei associated to H a (solid 
line) and P a (dashed line) during recombination. Among all tran- 
sitions listed in Table JTJ , only a few (one in case of H a , 3 in case 
of P a ) contribute to the final value of Ei . 



3 COSMOLOGICAL HYDROGEN 
RECOMBINATION 

The first scenario where resonant scattering may modify the 
pattern of the CMB polarisation anisotropies is precisely 
the cosmological hydrogen recombination. The CMB pho- 
tons observed today were essentially last scattered during 
the redshift range z £ [700,1300],, when the temperature 
of the Universe was su fficiently low to permit the forma- 
tion of neutr al atoms l|Zeldovich. Kurt fc Svunvaevl 1 19681 : 
|Peebledll968h . During the process of recombination of the 
hydrogen and helium atoms, there was a net generation 
of photons which introduced distortions to the CMB 
blackbody spectrum. Except for the far Wien tail of the 
CM B, these distortions are e xpected to be very small 
(e.g. iRubino-Martm et alj [200^ 1. with typical values of the 
order of AI/I w 10 -7 , making their detectio n chall enging. 
However, as studied in IRubino-Martm et alj l|2005l ). these 
lines leave also an imprint in the temperature angular 
power spectrum of the CMB. When observing the CMB at 
a frequency in the vicinity of a redshifted hydrogen line, the 
small optical depth of this transition influences the angular 
distribution of the CMB, due to the resonant scattering of 
the photons in this line. Therefore, we can use the optical 
depths connected with these transitions to explore the 
consequences of the overpopulation of the hydrogen atom 
levels at these epochs. 



Here, we build upon the work of IRubino-Martm et alj 

(2005) and compute the change in the CMB polarisation 
anisotropies during hydrogen recombination. In particular 
we focus on the main (An = 1) transitions of the Balmer 
and Paschen series, i.e. H a and P a : the first line pro- 
vides the highest amplitude, whereas for the second line 
the corresponding redshifted frequency lies in the frequency 
band where most of the future CMB experiments will ob- 
serve, and the foreground contamination is minimal. Note 
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Figure 5. Change in the TT (a panel), EE (b panel) and TE (c panel) power spectra due to resonant scattering on H a versus observing 
frequency. Solid, dashed and dotted lines correspond to multipoles I = 783, 890 and 1001, respectively. 




Figure 6. Same as in Figure J5}, but for the P a line. Note the different amplitude and frequency ranges. 



th at in addition to the b ound- bound transitions described 
in iRubino-Martfn et al.l (|2005l ). there is also a contribu- 
tion from the free-bound emission during recombination 
|Chluba fc Sunvaevl 1200(3 ). However, in the considered fre- 
quency band, their contribution to the total optical depth in 
the redshift range w here recombination occurs is very small 
l|Chluba et alj|2006h . 

For the H a line, there are seven permitted (electric 
dipolar) transitions from the third shell (n — 3) to the sec- 
ond one (n = 2), depending on the particular configuration 
of the electronic orbital, intrinsic (spin) and total angular 
momenta of the levels involved. For the P a case, the num- 
ber of transitions amounts to seventeen. In Table [T] we list 
these transitions. For each line, the total optical depth ob- 
served at the corresponding frequency will be the sum of all 
contributing transitions. For each of trans ition j — > i, the 
optical depth Tij is given by l|Sobolevlll94r3 ) 



where Aji is the sponteaneous emission Einstein coefficient, 
Xij is the wavelenght, and n%— j, gi-j are the populations 
and degeneracy factors of each of the two levels involved in 
the transition, respectively. H(z) is the Hubble function at 
redshift z. In our computations, the populations for the dif- 
ferent sublevels having different orbital angular momenta 
(2s, 2p, 3s, 3p, 3d, 4s, 4p, 4d and 4/) were tracked indepen- 
dently. Within each of these sublevels, we assummed that 
states having different total angular momenta j were evenly 
populated, i.e., as it is dictated by the Clebsch-Gordan coef- 
ficients and/or the statistical equilibrium condition. This is a 
reasonable assumption given that the energy differences be- 
tween the different fine-structure sub-states of a given shell 
are four order of magnitude smaller that its energy. Given 
that the recombination lines are formed when the number 
of photons (per ba ryon) which are able to exci te them falls 
below one (see e.g. iRubino-Martm et al.l (|2006l )), it is clear 
that the fine-structure sub-levels will maintain equilibrium 
due to the ambient CMB radiation field. Once we derive the 
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optical depth for each transition, we compute the average 
polarisation coefficient Ei associated to both H a and P a as 



E x = 



(28) 



with E\ 3 the polarisation coefficient for the particular 
transition i —> j. Figure Q shows the behavior of both 
Ei parameters through the recombination epoch in these 
circumstances: they show a somewhat constant value 
around ~ 0.35 for H a and ~ 0.27 for P a . This is mostly 
due to the presence of few dominant transitions (one for 
H a , three for P a ) having the largest weight in eq. (|28[) . In 
the case of H a , it is the P3/2 — > D5/2 transition driving the 
final value of Ei (~ 60% of the total r), whereas for P a 
these are D 7/2 -> F 7/2 (~ 40%), D 3/2 -> F 5/2 (~ 25%) and 
D 5/2 ^P 3/2 (~ 7%). 

Figure (J3J) displays the changes in the TT, EE and 
TE CMB angular power spectra versus observing frequency 
in the range where the H a is redshifted. The (remarkably 
smaller) amplitudes for the P a line are shown in figure ©. 
We have taken for the effective relative width of the line 
the value 5 x 10 -3 , (ax/vx, see eq. (|26[) ). Solid, dashed and 
dotted lines correspond to multipoles I = 783, 890 and 1001 
(corresponding for the maximum change in TE, TT and EE 
power spectra, respectively). These effects are considerably 
small, but however might be within the sensitivity limit of 
upcoming CMB experiments attempting to measure the B 
mode of polarisation anisotropies. Actually, when looking at 
H a , the amplitudes for the change in the TT and TE power 
spectra are comparable with the expected amplitude of the 
B modes at intermediate angular scales (Z 6 [60,600]) for 
a tensor to scalar ratio r of 0.03. Instead, the amplitudes 
for the changes in the EE power spectrum are closer to the 
same prediction for the B mode at larger angular scales, 
(I £ [10,30]). Also, the definite I and frequency dependence 
of these changes in the power spectra should make their 
detection more tractable in practice, although at those fre- 
quencies dust is suppossed to be dominating. In any case, 
we must conclude that, together with the characterisation 
of the B modes of CMB polarisation anisotropies, the detec- 
tion of such small signals remains challenging for ongoing 
and upcoming observational efforts. 




Figure 7. Relative fraction of ions for two different reionisation 
scenarios: in case (A) Pop III stars at 2 ~ 10 are able to ionise 
three times C, N and O only temporarily, so that those ions are 
absent until a new starburst at around 2 ~ 1 generates them 
again. This scenario is depicted by thin solid lines. In case (B), 
Pop III stars are able to keep a significant fraction of C, N and 
O ionised upto 2 ~ 1. In both cases, below 2 ~ 1, the star forma- 
tion rate drops considerably, and so do the abundances of highly 
ionised atoms. 




4 REIONISATION 

In this Section, we revisit the imprint of resonant scattering 
of CMB photons off metals produced at the end of the dark 
ages and the beginning of reionisation. We shall attempt 
to build a somewhat realistic model for the production of 
heavy elements during cosmic times, and to obtain from it 
an amplitude for the effect of scattering on metals produced 
by the first stars. 



4.1 Resonant Transitions 

A search for resonant transitions falling in a frequency range 
that permits probing the resonant interaction with CMB 
photons in the redshift range 2 g [5, 50] yields the list shown 



Figure 8. Volume average metallicity measured by a CMB ex- 
periment of angular resolution close to the arc-minute level and 
Au/u ~ 0.1. 



in Table[2] These are magnetic dipolar (Ml) forbidden tran- 
sitions, associated to elements and ions that should be more 
abundant during the epoch of cosmic dawn. After assum- 
ing a (volume-average) metal abundance of one part in one 
hundred compared to the solar value, we obtain typical val- 
ues for the optical depth in the range tx £ [io- 4 ,io~ 6 ]. 
I.e., we shall be working in the optically thin limit. Note 
that for the OI 63.2 line the polarisation coefficient is 
small (Ei — 0.01), since it will have its repercusion in the 
forthcoming discussion. 
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Table 2. Resonant transitions in atoms/ions expected to be produced during Reionisation 
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4.2 A Toy Model for the Enrichment History of 
the Universe 

In a similar way to BHMS, we shall consider two different 
models of reionisation: in case A we consider a reionisation 
scenario starting at z ~ 12, that is able to ionise twice 
a significant fraction of atoms like N, C and O only 
temporarily (down to redshift, say, z ~ 7). Only after a 
new starburst takes place in galaxies at z ~ 1 these ions 
appear again in the IGM. In case B, however, we observe 
the scenario where double ionised species are present in the 
IGM from z ~ 10 down to z ~ 0.5. The relative abundances 
of ions for C, N and O for this toy model are given in Figure 
0: thin lines correspond to case A scenario, thick ones to 
case B. These elements correspond to the most abundance 
species that are able to interact with CMB photons in the 
frequency range of interest for us. Neutral Carbon and 
Hydrogen dissappear as soon as first stars enter into scene 
for both reionisation models, but singly ionised species (CII, 
Oil, Nil) will eventually re-appear only if the flux of UV 
photons becomes lower than some threshold, (case A). The 
pattern of the polarisation anisotropics will be particularly 
sensitive to the presence or absence of certain species, as 
we next show. 



In our toy model, we shall assume that metals are 
located in the interior of haloes more massive than a given 
threshold mass, here taken to be M t h ~ 1O 8 M0. Inside 
these haloes the metallicites are taken to reach the level of 
5 times the solar metallicity. In practice, due to the finite 



angular and spectral resolution of the CMB experiment, 
future observations will be sensitive to the total amount 
of atoms/ions within a given cosmological volume probed 
by the CMB experiment, i.e., they will be sensitive to the 
volume average metallicity of the IGM. In Figure (JHJ we 
show the volume average metallicity to be measured by 
an experiment whose PSF and frequency response probes 
a region comparable to the Hubble radius at each epoch, 
as it should be the case for experiments having a spectral 
resolution of Av/u ~ 0.1. With this high dillution of the 
metallicities, the resulting amplitudes for the resonant 
scattering will be accordingly small. For simplicity, we 
shall restrict our analyses to the first ten lines considered 
in Table @, since they correspond to the strongest line 
of the most abundant elements. The amplitudes and the 
overall behaviour of the IGM metallicity given in Figure 
$Q are n °t t oo dif fe rent from obse r vation al estimates of 
ISchaye et all (l2003h: lAguirre et al.l (120041 ); iLedoux et ail 



|2003t ) or lProchaska et all (120031 ) 



The fact that all existing lines will introduce some 
anisotropies at redshifts dictated by the observing and res- 
onant frequencies translates into a non-trivial relation be- 
tween observing frequency and effective redshift of scatter- 
ing. In Figure (|U.) we display this effective scattering red- 
shift defined as 



(29) 



where the sum is over the resonant lines under considera- 
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Figure 9. (a) Effective redshift versus observing frequency, (b) Total optical depth versus redshift (c) Total optical depth versus observing 
frequency, (d) Average polarisation coefficient versus observing frequency. In all panels, thin solid lines refer to case A reionisation scenario, 
thick solid lines to case B. 



tion, and the optical depth for each transition is considered 
to be a function of the redshift corresponding to the observ- 
ing frequency. Again, the thin solid line corresponds to case 
A, thick solid line to case B. Although there is a general 
trend to have smaller effective redshifts with increasing ob- 
serving frequencies, the actual values may depend critically 
on the dominant lines at each frequency. This is the cause 
why panel (b) in Figure © (total optical depth Trot versus 
redshift) differs so much from panel (c), (tt ± versus ob- 
serving frequency): in panel (b) the distinction between the 
case A and case B scenarios looks much more feasible than 
in panel (c). As shown in panel (a), the case B scenario is 
observing, on average, at higher effective redshifts than case 
A, and this results in lower abundances and lower optical 
depths, making it more similar to the case A scenario. Both 
scenarios, case A and case B, could however be more easily 
distinguished if measurements of the effective polarisation 
coefficient were available; this coefficient, denned by 



p , s _ J2x T X [zx(v)] 

h i effW = — r — TTi — ■ 



(30) 



is plotted in Figure (]5Jl). The case B scenario (thick line) 
shows a profound dip at v ~ 140 GHz, and its origin 
is displayed in Figure (|10|l . This last figure provides the 
relative weight of each transition to the total optical depth 
at a given observing frequency: dotted lines correspond 
to Carbon, dashed lines to Nitrogen and solid lines to 
Oxygen. Increasing thickness corresponds to increasing 
line of entry on Table [5] for a given element, from top to 
bottom. We see that neutral Carbon (thin dotted lines) 
dominate the total optical depth at low frequencies, i.e., 
high redshifts. On the contrary, at low redshifts (or high 
frequencies) the optical depth is dominated by ionised 
Oxygen (thick solid line). At intermediate redshifts, the 
lines from Nil and CII compete (medium thickness dotted 
and dashed lines, respectively). The 63.2 /im 01 line (thin 
solid line) becomes prominent at v ~ 140 GHz, and in 
case B scenario where there is little Nil (since most of 
Nitrogen is in form of NIII) clearly dominates the total 
optical depth in a narrow redshift range. Since the polari- 
sation coefficient associated to this transition is relatively 
low (Ei = 0.01), it leaves a distinctive signature in the 
E 1 e g. Of course, these conclusions are consequences of our 



particular choice of reionisation models, and are not generic. 

Nevertheless, from a general point of view it is possible 
to state that the OI 63.2 fim transition will always leave 
the signature in the form of very low polarisation param- 
eter E 1 e g-, enabling its isolation from the rest of resonant 
lines. Further, it is a prediction for all models that the polar- 
isation coefficient should increase with increasing frequency, 
since in this regime it would be dominated by the ionised 
species OIII, Nil and CII. The measurement of polarisation 
anisotropies generated during reionisation via CMB scatter- 
ing off metals would not only provide a consistency check for 
the detection in the temperature/intensity side, but would 
also shed more light on how reionisation took place. How- 
ever, the detection of these polarisation anisotropies would 
require sensitivities of the order of 1-100 nK, and currently 
this is beyond the scope of most CMB polarisation projects. 



5 DISCUSSION AND CONCLUSIONS 

The measurement and interpretation of the polarisation 
intrinsic anisotropies of the CMB has remarkably tight- 
ened the cons t raints on several cosmological parameters, 
jSpergel et al.l |2006h . However, their characterisation is 
hampered by their low amplitude and the presence of 
contaminating signals. To date only anisotropies of the E 
mode of polarisation have been detected, but the detection 
of the B mode is currently the main goal of many ambitious 
CMB projects. The main difficulty that these projects will 
have to overcome is not only the extremely demanding level 
of sensitivity required to measure such a small signal, but 
also the exquisite control of systematics and foregrounds 
contaminants, the latter being usually orders of magnitude 
above the pursued component. 

In this context, we have presented a study of the 
secondary anisotropies of the CMB polarisation field caused 
via resonant scattering in two different scenarios: the cos- 
mological recombination and the cosmological reionisation. 
In both cases, the new polarisation anisotropies are E type, 
but of smaller amplitude than the intrinsic E type CMB 
polarisation fluctuations. In the recombination scenario, 
we have focused our analyses on the H a line, giving rise 
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Figure 10. Relative contribution of each line to the total opti- 
cal depth versus observing frequency. Dotted, dashed and solid 
lines correspond to Carbon, Nitrogen and Oxygen, respectively. 
Increasing line thickness corresponds to increasing line of entry 
in Table [2] 



to the largest signal among all other transitions. We have 
found a definite and distinctive pattern in multipole and 
frequency space for the change of the TE and EE power 
spectra. This effect takes place in the frequency range 
v G [350, 600] GHz, and is maximum at multipoles I ~ 900. 
The amplitudes of these signals are in general comparable 
to those of the intrinsic B-mode polarisation anisotropics. 
In a scenario where CMB experiments have the sensitivity 
and systematics control level required to characterise the 
B-mode anisotropics, the access to the study of the fluc- 
tuations introduced by the H a line should be then reachable. 

The polarisation anisotropies introduced by resonant 
scattering off metals during reionisation would provide a 
concluding consistency check for the measurements carried 
out in the intensity/temperature side. At the same time, 
they would shed some light on which transitions are respon- 
sible for most of the measured effect. In particular, the OI 
63.2 /xm line tends to lower the average polarisation coef- 
ficient Ei generated by all species, since it conserves very 
little memory of the polarisation of the scattered photon 
(£i,63.2nm = 0.01). On the other hand, the dominant tran- 
sitions associated to highly ionised species, like CIII and 
OIII, show relatively high values of E x , (Ei £ [0.5, 1]). The 
requirements for detecting these signals are, however, even 
more demanding than in the recombination scenario: typ- 
ical amplitudes are closer to the nK level than to the fiK 



level, even below the amplitudes expected for the primordial 
B-mode polarisation anisotropies. The final detectability of 
this signal will be critically associated with the detectability 
of the primordial B-mode anisotropies: if the latter, at the 
end of the day, can be characterised with some given accu- 
racy, then the possibility of pursuing signals that are still 5 
to 50 times smaller could be formally addressed. 
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APPENDIX A: THE DRAG FORCE 

A net flow of CMB photons will exert some drag force on 
a fluid of atoms that are able to interact with the photon 
gas via some resonant transition X, with drag force defined 
as " rate of momentum exchange per unit time and unit vol- 
umen" . Here we detail the calculation of the drag force ex- 
erted on the atom fluid for a given frequency of the CMB 
field. The expression will be integrated in atom momentum 
and photon frequency to yield the well-known result that 
the drag force will be proportional to fo(vx, rD^iVxnxirj), 
with nx (rj) the average comoving number density of atoms 
taking part in the resonance. 

For a given transition X, the optical depth for scattering 
at a given frequency is given by 



*(„) = 



9-u 



8nv x gi 



A x $(v) 



l + £«*P(-& 



(Al) 



where Ax is the Einstein coefficient for spontaneous emis- 
sion of the transition, g u and gi the statistical weights for 
the upper and lower levels of the transition, respectively, and 
T ex represents the excitation temperature, which in case of 
equilibrium with the radiation field would correspond to the 
CMB temperature. Finally, &(v) represents the line profile. 

The photon flow and drag force are taken to be paral- 
lel to the polar axis in the integration on q, (i.e., n). Only 
the dipole term fo(q, r;)^'i(k, q, rj) will contribute in the mo- 
mentum integration, and hence the drag force amplitude, 
expressed in terms of unit volumen occupied by the atom 
fluid, can be written as 



F 



(-Am) / dqq 2 f o (q,v)^i0<-,q,v) Q ca x 



dp fx(p,v)- 
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For simplicity, we have introduced a frequency-independent 
cross-section ax, which in practice is equivalent to approx- 
imate the profile function by a top-hat function of a cer- 
tain width Av. Thus, the expression for ax is given by 
a(v)/(Av$(v)). Note that / *i = Yf = ifoM corresponds 
to the sum over polarisation states of the first momentum 
of the photon distribution function expansion introduced 
in eq.(|3}. The atom fluid is assumed to be described by a 
Maxwell-Boltzmann distribution function (fx) in thermal 
equilibrium with the CMB photon field: 



fx(p,v) = 



n x (ri) 



[2nm x k B T} 



3/2 



exp ( 



~p 2 /[2nmxk B T]), (A3) 



with T — To(l + z\rj\) the CMB temperature monopole. 
The dominium of integration E is defined by the subset of 
possible values of p and /i (dp = p 2 dpdfid(p) verifying the 
resonance condition 



p ■ n = mxc 
that is, 



(A4) 



p n G 



. , , , vx .A x 
mm I mxc ( 1 ± 

v 4-7T 



, , . vx , A x 
max mxc 1 ± 

v 4-7T 



(A5) 



with Ax /2-7T the intrinsic natural width of the line. Since the 
integrand does not explicitely depend on fj,, we can integrate 
on this variable, replacing dfi by Afi — mxc Ax / (2nvp). 
Note that A^t ^ 2, and so 



m x c A x 



(A6) 



1 2 2tvv 

Note, as well, that from the resonance condition of eq, ()A4[l 
we must have 



p > mxc(l--^], (v>v x ), 



p > m x c[ 1 

v 



(v < vx). 



(A7) 
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Therefore, the minimum value of p allowed in the integra- 
tion is 

m x c A x 



PminM = max 



2 2w 



± mxc[ 1 — 



vx 



(A9) 



The first value will be used only if | vx /v — 1 ^ Ax /(&irv x ), 
and will give a negligible contribution, as we shall see. The 
maximum value for p will be pmax = mxc. However, it will 
be approximated by infinity in the integral, since mxc 2> 



\/2raxkBT in our case. After introducing Afi in eq. (IA2[) . 
trivially integrating in d<f>, and integrating in p one obtains 
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A, 



mxc 



2ttv ^2m x k B T 

exp(-p 2 min (v)/(2mk B T)). (A10) 

This equation shows that, for a given frequency, the fraction 
of atoms being affected by the drag of CMB photons is 



rx = 



m x c 



2tvv ^/2mxk B T 



^p(^p 2 min (v)/(2mk B T)). (All) 



^,From eq. (|A9|l . it is clear that the largest contribution of 
the integral in eq. (|A10[l comes from the vicinity of q — qx = 
hvx /c, and this allows rewriting most of the terms out of 
the q integral, yielding 



(-4m)q% (/o*i)sx Qx (c a x n x ) 



Ay, 



mxc 



2ttvx ^2m x k B T 



line 



dq exp 
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We next sweep the vicinity of the line, and split this inte- 
gral in three separate ones, according to the three different 
possible values of p m { n (v). The final result is: 



(-4-Ki)(ca x n x ) ~^:<& (/o*i)« 



qx 



erfc 



(—) 



Ay 



exp - 



mxc A x V 
^/2rnxk B T ^vx J 



(A13) 



The first term in brackets accounts for the integral through- 
out the thermal width of the line. The argument of the com- 
plementary error function is very close to zero, so this term 
is very close to unity. That is, practically all atoms are ex- 
periencing some drag force exerted by the CMB photons. 
The second term in brackets provides the integral within 
the natural width of the line, and its contribution can be 
safely neglected. Hence 

~ (-47rt)(c axn x ) ^q x (/o*i)« qx . (A14) 

I.e., the momentum qx = hvx/c is transferred to all atoms 
at a rate given by c axnx in a momentum shell of radius 
qx and thickness hAx / (2-kc) . 



